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SECTION  I 
INTRODUCTION 

The  performance  of  airborne  electrooptical  systems  may  be  seriously 
affected  by  the  angular  as  well  as  the  linear  vibration.  Unlike  the 
linear  vibration,  very  little  experimental  data  are  available  for  the 
angular  vibration.  The  purpose  of  this  report  is  to  investigate  the 
relationship  between  the  linear  and  angular  structural  vibration 
characteristics  thereby  developing  prediction  schemes  for  the  angular 
vibration  amplitude  from  a prescribed  l.inear  vibration  environment. 
Precisely  speaking,  the  linear  vibration  refers  to  the  transverse 
deflection,  whereas  the  angular  vibration  is  due  to  the  rotary  mode  of 
internal  bending,  the  effect  of  which  is  included  in  the  Timoshenko 
beam  equation  (Reference  1).  For  a small  amplitude  vibration,  the  angular 
(rotary)  deflection  can  be  related  to  the  slope  of  the  linear  (transverse) 
displacement.  This  provides  a theoretical  basis  for  relating  the  linear 
and  angular  vibrations. 

As  a theoretical  model.  Section  II  considers  the  Bernoulli -Euler 
beam  with  simply  supported  ends  and  formulates  the  stochastic  dynamics 
for  the  linear  and  angular  deflections.  For  maximal  simplicity,  an 
idealized  random  forcing  function  is  imposed  to  generate  the  stationary 
response  dynamics.  Since  the  stochastic  beam  analysis  has  already 
been  presented  (References  2 and  3),  only  the  main  steps  of  analysis 
which  yield  the  linear  and  angular  variances  as  infinite  series  will  be 
summarized.  It  is  important  to  note  that  the  response  variance  has 
the  same  form  for  both  the  beam  and  plate  problems  when  the  normal  - 
mode  representation  is  used.  This  is  why  a beam  analysis  can  provide 
the  basic  relationship  of  stochastic  dynamics  which  is  also  valid  for 
the  plate  structure.  To  be  realistic,  the  aircraft  fuselage  must  be 
modeled  by  a stringer  stiffened  cylindrical  shell.  Since  such  a model 
is  not  amenable  to  qualitative  analysis,  an  alternate  approach  would 
be  to  build  up  the  aircraft  fuselage  by  piece-wise  flat  panels  joined 
together  by  stiffeners.  The  motivation  for  the  present  analysis  is  that 
a beam/plate  model  can  describe  certain  essential  features  of  aircraft 
fuselage  vibration  in  a qualitative  manner. 
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In  Section  III,  an  attempt  is  made  to  derive  a qualitative  relation- 
ship between  the  linear  and  angular  rms  amplitudes,  which  is  valid  for 
the  general  vibration  environment.  This  has  been  accomplished  by 
averaging  the  linear  and  angular  variances  over  the  beam,  thereby 
integrating  away  the  normal  modes  on  which  the  analysis  of  Section  II 
is  based.  By  truncating  the  variances  to  lowest  orders,  a simple 
prediction  (Equation  27)  has  been  obtained  to  relate  the  beam-averaged 
linear  and  angular  rms  amplitudes.  However,  due  to  the  scarcity  of 
angular  vibration  data,  comparisons  of  predicted  and  measured  angular 
vibration  are  necessarily  quite  limited. 

Because  of  the  difficulties  associated  with  many  angular  vibration 
measurement  techniques,  a practical  approach  is  to  calculate  a discrete 
angle  by  differencing  two  linear  accelerometer  signals  and  then  dividing 
by  the  separation  distance.  The  main  goal  of  Section  IV  is  to  provide 
a prediction  for  the  discrete  angular  vibration.  Based  on  some  laboratory 
test  data,  the  accuracy  of  the  prediction  scheme  as  given  by  Equation  33 
is  about  ± 20%  of  the  measured  angular  vibration.  Although  quite  crude, 
it  bolsters  one's  confidence  in  that  the  prediction  scheme  is  intended 
for  qualitative  purposes  and  the  present  vibration  data  are  limited 
only  to  a gross  confirmation  of  predictions. 

The  discrete  angle  defined  by  the  difference  of  two  linear  acceler- 
ometer signals  is  not  well  suited  for  statistical  analyses,  especially 
when  the  separation  distance  becomes  small.  This  is  because  the 
differencing  operation  tends  to  accentuate  the  noise  part  of  random 
signals.  Consequently,  the  smaller  the  separation  distance,  the  more 
acute  the  differencing  error.  In  Section  V an  alternate  definition 
is  presented  of  the  discrete  angular  vibration  which  involves  the 
difference  of  elements  of  the  covariance  matrix  and  not  the  random 
signals  themselves.  Surprisingly,  such  a definition  of  the  angular 
motion  has  not  yet  been  explored.  Its  practical  implication  is  being 
investigated  and  will  be  discussed  in  a later  report. 
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SECTION  II 


THEORETICAL  ANALYSIS 


The  transverse  deflection  of  a Bernoul 1 i -Euler  beam  (Figure  1)  is 
governed  by  the  equation 


where  El  is  the  flexural  rigidity,  m is  the  mass/unit  length  of  the 
prismatic  beam,  B is  the  damping  constant,  and  f(x,t)  is  the 
excitation/mass.  It  will  be  assumed  that  the  forcing  function  splits 
into  the  spatial  part  (q(x))  and  the  temporal  part  (g(t) ) 


As  pointed  out  in  the  Introduction,  the  Timoshenko  beam  incorporates  the 
angular  effect  of  internal  bending,  whereas  the  Bernoul 1 i-Euler  beam 
does  not.  Therefore,  it  might  appear  that  the  use  of  Equation  1 is  not 
appropriate  for  angular  vibration  work.  As  shown  in  the  Appendix, 
however,  both  beam  theories  give  the  same  prediction  scheme;  hence  it  is 
more  economical  to  base  the  analysis  on  Equation  1,  which  is  somewhat 
simpler  than  the  Timoshenko  beam  equation.  The  stochastic  dynamics 
of  Equation  I have  already  been  investigated  by  Eringen  (Reference  2), 
and  Bogdanoff  and  Goldberg  (Reference  3),  so  that  it  suffices  to 
outline  only  the  main  steps. 


(1) 


f (x,t)  = q(x)  g (t). 


Y ( X, t ) 


X 


Figure  1.  Simply  Supported  Beam 
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1.  LINEAR  AND  ANGULAR  DEFLECTIONS 

For  the  simply  supported  beam  of  length  L,  the  undamped  normal 
mode  sin(mrx/L)  can  be  used  to  express  the  spatial  behavior 

y(x»t)=Z  yn  (t)  sin  ( JlZLi  \ (2a) 

n = l L 


q(x)  = I qn  sin  (^) 
n = I 


(2b) 


The  insertion  of  Equation  2 into  Equation  1 gives 

yn  +2/3  yn  +o>2n  yn  = qn  g(t), 

? 4 

where  the  dot  denotes  3/3t,  and  = (nTr/L)  El/m.  The  particular 
solution  with  zero  initial  conditions  is 

r\ 

yn(t)  = dn  Jtohn(,-T)9(r)dr. 

Flere  the  impulse-response  function  is  given  by 

exp(-/3t)  , , 

hn  (t)  = — sin  a>n  t 


(3) 


(4) 


where  w 


7f~Z  2 

— 3 . By  letting  t^-*-  - one  may  rewrite  Equation  3 as 

yn(')=qnf\lr)  g(t-r)dr  (5) 


The  introduction  of  Equation  5 into  Equation  2a  would  then  yield  the 
linear  deflection 


y(x,t) = £ 
n = l 


sin 


( !!ZpL)  qn  J hn(r)  g(t-r)  dr, 


(6) 


subject  to  the  arbitrary  forcing.  For  a small  amplitude  vibration,  the 
angular  deflection  0(x,t)  is  approximated  by  the  derivative  dy(x,t)/dx, 
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Thus,  the  linear  and  angular  deflections  are  formally  expressed  by 
Equations  6 and  7 respectively. 

2.  VIBRATION  AMPLITUDE 

When  the  forcing  function  is  random,  so  are  the  linear  and  angular 
deflections.  Let  us  therefore  compute  the  variance  of  y(x,t) 


where  < > denotes  statistical  average  or  expectation.  It  must  be 

pointed  out  that  the  statistical  independence  of  q(x)  and  g(t)  has  been 
invoked  in  Equation  8,  which  greatly  simplifies  the  subsequent  analysis. 
For  instance,  this  assumption  cannot  be  applied  when  a beam  or  plate 
is  excited  by  convecting  turbulent  pressure  fluctuations  (References  4 
and  5).  For  stationary  random  forcing,  the  correlation  < q(t-r^ JgCt-x^)  > 
is  a function  of  the  time  difference  x = - x-j  . By  the  Wiener-Khinchin 

theorem  (Reference  6),  the  stationary  correlation  of  the  forcing 
function  denoted  by  R (x)  is  related  to  the  power  spectral  density  (psd) 
denoted  by  Fg(co). 


00 


hn  (T|)  hk  (t2)  <g  (t  - T|)  g(t  - r2)>  dr,  dr2, 


(8) 


0 0 


oo 


- oo 


(9) 
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Under  the  stationary  forcing,  it  is  evident  that  Equation  8 is  inde- 
pendent of  t,  hence 

n,k 


(10) 


where 


CD  GO  CD 

lnk  = i / do.  Fg(co)  / h (r,)  e-icur,  dr,  J hk  (r2) e 'WTz  dr2  . 

0 0 i 1 1 


Similarly,  the  stationary  angular  variance  is  given  by 

<02(*)>=  (^)  I nkcos(^)cos(~)<qnqk>  lnk 
n,k 

With  the  use  of  Equation  4,  one  finds  by  a simple  computation  that 


(12) 


/ hn(r,)  exp(-icuT,)  dr,  = 
J0 


cut,  + 2\/3<jj-uj‘ 


= Hn  (a>) , 


which  is  the  frequency-response  function.  Hence,  Equation  11  may  be 
put  in  the  compact  form 


oo 

I nk  = I J Fg  ( cjli)  Hn  ( cu)  Hk(cu)  dcu. 


(13) 


where  * denotes  the  complex  conjugate.  When  the  normal -mode  repre- 
sentation is  used  (Reference  2),  the  form  of  linear  and  angular  variances 
is  the  same  for  both  the  beam  and  plate  problems  in  that  it  consists 
of  the  normal  modes,  the  spatial  correlation  < qnqk  >,  and  the  1^. 


To  explicitly  evaluate  Equation  13,  one  must  further  specify  Fg(w). 
For  a typical  R^t)  of  the  exponential  form,  the  corresponding  psd  has 
the  form 


Fg  (OJ)  = 


^0^0 
U)q  +OI2 


» 


(14) 
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where  Fq  is  a scaling  factor,  and  coq  is  the  relaxation  constant.  In  the 
limit  as  Wq  °°,  Equation  14  reduces  to  the  white-noise  spectrum 

Fg  (cu)  = F0  . (15) 

For  simplicity,  let  us  evaluate  Equation  13  under  the  white-noise 
forcing  Equation  15.  The  result  by  contour  integration  is 

4tt£f0 


nk  ~ 


(o>n  - “4/32)2+(4/3w'n)‘ 


(16) 


which  is  2iri  times  the  residues  at  co  = i$  ± w'.  For  n = k,  it  simplifies 
to 


ttF0 

1 nn  “ „ a z 


4/9  0*5 

4 

From  the  definition  of  cjn,  it  is  seen  that  Inn  ^ 1/n  , hence  the 

magnitudes  of  I^n  fall  off  rapidly  with  increasing  n.  Magnitudes  of 

I , (n/k)  with  respect  to  I will  be  assessed  in  Section  1 1 -4 . 
nk  v ' K nn 


(17) 


3.  RANDOM  FORCING  q(x) 


In  view  of  Equation  2b,  the  mean  value  < q^^  > has  the  form 

«,'*.»'-»(*?)-"(¥)•  (18) 

G O 

where  R (x,y)  = < q(x)q(y)  > is  the  spatial  forcing  correlation.  When 
spatial  homogeneity  is  assumed,  Rq(x,y)  is  a function  of  the  separation 
distance  only.  By  letting  z = x - y,  Equation  18  can  then  be  put  in  the 
form 

= T f0L  Rq(z)  sin(nFL~')dz 

when  R (z)  is  assumed  fairly  sharp  (Reference  7).  Now  expand  sin(mr(z+y)/ 
L)  = sin(mrz/L)  cos(mry/L)  + cos(mrz/L)  sin(nTry/L).  Since  Rq(z)  is  even, 
it  is  seen  that  Rq(z)  sin(mrz/L)  integrates  out  to  zero,  and  hence 
Equation  18  becomes 


< Qn  ^k>=  ^nk  f Rq  ( z)  cos  ( L ) dz» 


(19) 
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where  the  Kronecker  6pk  =1  for  n = k and  = 0 for  n t k.  For  the 
infinitely  sharp  R (z)  = K 6(z)  (K  being  a constant).  Equation  19 
reduces  to 

<qnqk>  = KSnk  (20) 

Under  the  random  q(x)  with  no  spatial  correlation,  the  variances  shown 
in  Equations  10  and  12  therefore  become 

<y2(x)>=  Z sin2(^p)Klnn,  (21a) 

n = l 


< 02  (x)  > 


(21b) 


-4 

Because  I ~ n , the  series  representations  are  convergent,  even 
though  the  total  power  of  f(x,t)  is  infinite.  To  be  more  physically 
realistic,  the  forcing  correlation  must  exhibit  both  the  spatial  and 
temporal  relaxations.  Although  this  can  easily  be  incorporated  by 
Equations  14  and  19,  such  is  not  needed  for  the  present  purpose. 


4.  LOCALIZED  FORCING  q(x) 

In  a laboratory  test  situation,  it  is  most  convenient  to  apply  the 
external  forcing  to  certain  points  on  a structure.  In  such  a case,  the 
< qnqk  > tenn  is  simply  a qnq ^ . When  Equation  20  is  not  invoked,  the 
variances  shown  in  Equations  10  and  12  are  double  sums  over  n and  k. 

For  a slightly  damped  beam  it  will  be  shown  that  I are  orders  of 
magnitude  larger  than  I^k  (n  f k);  hence,  the  variances  can  still  be 
reduced  to  a single  sum  even  when  q(x)  is  nonrandom.  To  this  end, 
some  lower-order  Ink  were  computed  from  Equation  16  with  6 = co^/10,  and 
their  relative  magnitudes  have  been  summarized  in  Table  1. 
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TABLE  1 

RELATIVE  MAGNITUDES  OF  THE  LOWER-ORDER  Ink 


Tnk 

!11 

l22 

l33 

I =1 
x12  *21 

^^l 

1 32^23 

Magnitude 

1 

1/16 

1/81 

1/1406 

1/40375 

1/26406 

From  the  table, 
slightly  damped 


one  may  therefore  conclude-  that  'nk  ~ n Snk 
structure  under  the  white  noise  excitation. 


for  a 
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SECTION  III 

THE  BEAM-AVERAGED  VIBRATION  RMS  AMPLITUDE 


Since  Inn~n'^,  the  first-term  approximation  of  Equation  21a 
introduces  an  error  on  the  order  of  (1/16) 


< y2(x)  > = sin2  KIU.  (22) 

2 

Because  of  the  factor  n , on  the  other  hand,  the  first-term  truncation  of 
Equation  21b 


< 02(x)  > = 


(23) 


discards  a term  of  the  order  of  (1/4),  which  is  not  acceptable.  To 
the  same  order  of  approximation  as  Equation  22,  the  second  and  third 
terms  must  then  be  included 


< 82l.)  > = (-J-)2  n|i  o2  COS2  (-ip.)  KI„„.  (24) 

2 2 

For  a localized  forcing  at  the  mid-span  (L/2),  K~q  ~ sin  (mr/2),  so 
that  the  second  term  (n  = 2)  will  then  be  absent  in  Equation  24. 


For  the  overall  description,  the  variances  will  now  be  averaged  over 
the  beam.  First,  by  applying  L dx  to  both  sides  of  Equations  22, 

Vrri  [ 2 

L / < y (x)  > dx,  has  the 

following  form 


rms 


V 


KIll'2 


(25) 


Next,  averaging  Equation  24  over  the  beam  in  a similar  manner,  one 

4 

obtains  with  the  use  of  I = I^/n 
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e 


rms" 


7 r 
L 


(26) 


/ 1 l 9 

where  6 = yL  / < 0 (x)  > dx  is  the  beam-averaged  rms  of 

rms  o 

angular  fluctuations.  Note  that  by  averaging  over  the  beam  we  have 
eliminated  the  explicit  dependence  of  Eguations  22  and  24  on  the  normal 
modes,  sin(mrx/L),  which  reflect  the  simply  supported  end  condition. 

It  is  therefore  anticipated  that  Equations  25  and  26  would  apply  to  any 
beam,  at  least  qualitatively,  regardless  of  the  end  conditions.  The 
practical  implication  of  beam  averaging  is  that  it  frees  the  formulation 
from  a particular  end  condition.  This  is  important  in  applications 
because  the  boundary  (or  end)  conditions  of  a structure  embedded  in 
airframes  cannot  be  specified  precisely.  By  combining  Equations  25  and 
26,  the  main  result  of  this  section  is 


where 


®rms  = A ( L ) yrms  ’ 


(27) 


A = 


( l l/n2)1/2  = 1.17, 
n=  I 


for  the  first-term  truncation 
Equation  23 

for  the  random  q(x). 


l (1  + 1/8)1/2  - 1.06, 


for  the  localized  q(x)  at  L/2. 


Note  that  only  the  lowest  order  modes  contribute  explicitly  to  the  above 
relation.  Although  Equation  27  was  obtained  from  the  beam  analysis,  it 
should  also  apply  to  a plate  along  the  coordinate  axes.  This  is  because 
the  basic  structure  of  variance  is  the  same  for  both  the  beam  and  plate 
problems  when  the  normal -mode  representation  is  used. 


Due  to  the  difficulty  in  angular  measurement  techniques,  angular 
vibration  data  are  very  scarce.  Also,  it  is  almost  impossible  to  find 
a matching  set  of  linear  and  angular  vibration  data  taken  concurrently  at 
the  same  location  on  an  airframe.  Some  angular  vibration  psd  have  been 
reported  by  Sher  (Reference  8)  during  the  flight  tests  of  an  NC-135 
aircraft.  Later,  Davis  (Reference  9)  measured  the  linear  vibration  psd 
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on  another  aircraft  of  the  same  type.  Apart  from  the  apparent  difference 

in  aircraft  and  flight  conditions,  it  is  most  unfortunate  that  the 

transducer  locations  of  Davis  are  not  the  same  as  those  of  Sher.  In 

spite  of  this,  an  attempt  has  been  made  to  select  a set  of  linear  and 

angular  psd  data  whose  transducer  locations  are  as  close  as  possible. 

For  the  present  comparison,  Figure  30  of  Davis  (Reference  9)  and  Figure 

8 of  Sher  (Reference  8)  were  chosen  as  a typical  set  of  the  linear  and 

angular  psd's.  They  are  reproduced  here  in  Figures  2 and  3,  respectively. 

To  estimate  the  rms  of  vertical  displacements,  one  may  approximate  the 

-5  -2.7 

psd  of  Figure  2 by  a power  formula  of  the  form,  (5x10  ) to  , which 

in  essence  replaces  the  psd  by  a straight  line  in  the  figure.  Then, 
analytically  integrating  this  power  formula  over  the  frequency  range 
of  [2,°°],  one  obtains  yrm$  - 3 x 10  in.  Therefore,  by  treating  the 
entire  fuselage  as  a beam-like  structure,  L = 2100  in..  Equation  27  with 
A = 1.17  yields  the  rms  of  angular  vibration  6rms  “ 5.24  x 10  ^ rad. 

Now,  to  provide  an  experimental  check,  0rms  can  be  estimated  directly 
from  Figure  3.  Again,  for  convenience,  the  psd  of  Figure  3 may  be 
approximated  by  a three-piece  power  formula  as  depicted  in  the  figure. 

Then  integrating  this  power  formula  over  the  appropriate  frequency  range, 
the  measured  rms  of  angular  vibration  is  estimated  to  be  0rms  - 6.23  x 10 
rad.  Although  the  measured  0rms  is  almost  20%  higher  than  the  predicted, 
such  a comparison  is  very  encouraging  in  that  linear  and  angular 
vibration  data  used  here  have  limited  correspondence.  The  scarcity  of 
vibration  data  impedes  further  confirmation  of  the  prediction  scheme.^ 
Therefore,  an  alternate  prediction  will  be  developed  in  the  next  section 
in  terms  of  the  discrete  angle,  which  is  readily  measurable  by  the 
conventional  accelerometers. 


1.  Some  additional  sets  of  linear  and  angular  vibration  data  have  been 
supplied  by  W.  R.  Davis,  Jr.  Although  these  data  seem  to  support 
the  prediction  capability  of  Equation  27  within  + 20%,  the  estimation 
of  rms  values  from  the  raw  test  data  involves  a number  of  gross 
approximations.  Hence,  the  details  of  comparison  will  not  be  presented 
here. 
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Figure  2.  PSD  of  Transverse  Vibration 
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CVJ 
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Figure  3.  PSD  of  Angular  Vibration 
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SECTION  IV 
THE  DISCRETE  ANGLE 


Conventional  transducer  technology  does  not  permit  a direct 
measurement  of  angular  vibration  over  a wide  range  of  frequencies.  In 
practice,  the  angle  is  approximated  by  differencing  two  linear  transducer 
signals  and  dividing  by  the  separation  distance.  That  is,  the  discrete 
angle  is  defined  by 

Ay  »(»!■«)  -»(«..»)  (28) 
Ax  Ax 

where  Ax  = - x^.  The  discrete  angle  will  approach  the  true  angle 

0(x,t)  as  Ax  ■+•  0.  In  the  actual  measurement,  Ax  is  anything  but 
infinitesimal,  hence  Equation  28  is  a working  definition  of  a measurable 
angle  rather  than  a limiting  process  to  the  true  angle. 

By  repeating  the  stochastic  dynamic  analysis  of  Section  II,  the 
variance  of  (Ay/Ax)  can  readily  be  computed  and  its  first-term  trunction 
is 


With  the  use  of  the  transformation  x^  = x +Ax/2  and  x^  = x -Ax/2,  one 
has  sin (ttx1  /L ) - sin(Trx2/L)  = 2 cos(itx/L)  sin(TrAx/2L) , hence 


Again,  averaging  Equation  30  over  the  beam,  the  beam-averaged  rms  of  the 

y_i  r 2 

L / <(Ay(x)/Ax)  >dx,  has  the  following 
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form  similar  to  Equations  25  and  26 


Since  (2  sin(TrAx/2L)/Ax)  = tt/L  as  Ax  -*  0,  it  is  seen  that  Equation  31 

3 2 1 /2 

reduces  to  Equation  26  with  the  factor  ( l 1/nc)  7 replaced  by  unity. 

n=l 

The  combination  of  Equations  25  and  31  would  then  yield  a relationship 
between  the  linear  and  discrete  angular  rms  of  vibration  amplitudes 


As  pointed  out  before,  the  factor  (2  sin(irAx/2L)/Ax)  has  the  value  (tt/L) 
in  the  limit  as  Ax  0.  In  the  other  extreme  case  of  Ax  = L,  however,  it 
has  the  value  (2/L).  Note  that  this  factor  varies  almost  57%  between  the 
two  extreme  cases.  In  practice,  the  separation  distance  is  neither 
infinitesimal  nor  the  entire  L;  let  us  therefore  choose  Ax  = L/2,  for 
which  (2  sin(7rAx/2L)/Ax)  has  the  value  (2v/2"/L).  In  terms  of  the  Ax, 
Equation  32  then  becomes 


This  is  a counterpart  of  Equation  27  for  the  discrete  angle  case.  The 
difference,  however,  is  that  Equation  33  refers  to  the  separation 
distance  Ax  and  not  the  beam  length  L.  Like  Equation  27,  the  relation 
shown  in  Equation  33  does  not  depend  on  the  particular  shape  of  normal 
modes;  the  details  of  spatial  characterization  have  been  eliminated  by 
the  beam  averaging.  In  this  respect.  Equation  33  is  an  overall 
description  of  beam  vibration  in  terms  of  the  primary  mode,  which  may 
qualitatively  be  valid  regardless  of  the  end  condition. 

Davis  and  Guckian  (Reference  10)  have  reported  some  measurement  of 
(Ay/Ax)  on  a 156"  x 60"  x 30"  bench  plate  externally  excited  to 
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simulate  flight  conditions.  They  have  presented  both  the  linear  and 
angular  vibration  data  in  terms  of  the  transfer-function  magnitude. 
Unfortunately,  they  did  not  report  the  linear  vibration  data  of  the 
accelerometers  from  which  the  angular  deflection  was  derived.  It  is 
therefore  necessary  to  pick  out  linear  vibration  data  from  their  report 
in  such  a way  that  the  linear  and  angular  vibration  data  would  have  the 
closest  correspondence.  For  the  angular  vibration  data  given  in 
Figures  54  through  57  of  their  report,  the  linear  transfer-function 
magnitudes  of  their  Figures  33  and  35  were  chosen  as  the  matching  linear 
vibration  data.  To  estimate  the  rms  of  fluctuations  from  a transfer- 
function  magnitude,  the  following  procedure  has  been  adopted;  (i)  Assume 
a forcing  psd  of  the  form  Cw  2 (C  being  a constant),  (ii)  approximate  the 
transfer-function  magnitude  by  a piecewise  power  formula,  similar  to 
approximating  the  psd  by  a power  formula  in  Section  III,  (iii)  multiply 
the  squared  transfer-function  magnitude  by  the  forcing  psd,  and  (iv) 
integrate  over  the  appropriate  frequency  range.  For  each  of  the  angular 
vibration  data,  the  measured  (Ay/Ax)^ms  is  estimated  by  the  above  procedure 
and  compared  with  the  predicted  value  based  on  Equation  33,  as  summarized 
in  Table  2. 

Discounting  the  case  of  Roll  #2  whose  measurement  is  apparently 
contaminated,  one  can  predict  by  Equation  33  the  measured  angular  rms 
within  less  than  + 20%.  Although  this  is  a fairly  wide  latitude,  it 
must  be  remembered  that  the  comparison  of  Table  2 is  subject  to  several 
limitations.  First  of  all,  the  prediction  scheme  (Equation  33)  is 
valid  only  qualitatively;  that  it  can  work  so  satisfactorily  is  indeed 
a pleasant  surprise.  Second,  the  sets  of  linear  and  angular  vibration 
data  do  not  match  exactly  with  respect  to  the  transducer  locations. 

Finally,  the  values  of  yms  and  (Ay/Ax)rm$  estimated  from  the  transfer- 
function  magnitude  data  are  only  approximate. 
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SECTION  V 

PRACTICAL  DEFINITION  OF  THE  DISCRETE  ANGULAR  VIBRATION 

? 

According  to  Equation  28,  the  direct  measurement  of  < (Ay/Ax)  > 

involves  taking  the  difference  of  the  signals  of  two  accelerometers  which 

are  separated  by  a predetermined  distance  Ax.  In  practice,  the 

differencing  of  random  signals  is  not  a reliable  operation  because  it 

tends  to  accentuate  the  noise  part  of  random  signals.  Clearly,  the 

smaller  the  separation  distance  the  more  acute  the  differencing  error. 

An  alternate  definition  will  therefore  be  presented  of  the  discrete 

angular  vibration  amplitude  which  involves  subtraction  of  the  well-behaved 

statistical  averages  and  not  the  random  signals  themselves.  To  this  end, 

square  the  right-hand  side  of  Equation  28  and  carry  out  the  term-by-term 

averaging.  The  discrete  angular  variance  is  then  expressed  in  terms  of 

the  linear  variance <y  (x.)  > and  cross-correlation  < y(x,)  y (x0)  > 

/ . \ 2 1 1 ^ 

Ay  (*i , xA) 


Ax 


> = 


J^Pjz  [<v2(«,)>  + <r2(*2)>-2 

<»(*,)y(«2)>]- 


(34) 


Note  that  < y (x.j)  > and  < y(x^)  y (x^)  > are  respectively  the  diagonal 

and  off-diagonal  elements  of  the  covariance  matrix  < y(x.. ) y (x^.)  > (i,  j 

= 1,2).  Hence,  in  terms  of  the  covariance  matrix.  Equation  34  is  a 

2 

practical  prescription  for  computing  < (Ay/Ax)  >,  especially,  when  Ax  is 
smal 1 . 


Although  the  variance  was  computed  in  Section  II,  the  cross-correlation 
has  not  yet  been  discussed.  For  the  simply  supported  Bernoul 1 i-Euler 
beam,  < yU-jJy^)  > can  be  written  down  immediately  by  replacing  one  of 
the  y(x)'s  in  Equation  21a  by  y(x-j)  and  the  other  y(x)  by  y^) 


y (*i)y  (*2) 


I 

n = I 


sin 


(n7Tx,  \ / n7Tx2  \ 

— L— ) S'n  \ L ) 


K I 


nn 


(35) 
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With  the  use  of  x]  = x + Ax/2  and  x2  = x - Ax/2,  Equation  35  becomes 


<»(«,)»('!)>■  £ - c»s2  (-“f*-)]*1™-  (36) 

n = I x 

As  a consistency  check,  the  insertion  of  Equation  36  together  with  the 

corresponding  < y (x^)  > into  Equation  34  reproduces  the  previously 

defined  angular  vibration  amplitude,  the  leading  term  of  which  is 

Equation  30.  This  shows  that  Equation  34  is  just  another  definition  in 

terms  of  the  covariance  matrix.  Recall  that  the  assumption  of  the 

infinitely  sharp  forcing  correlation,  < q(x^)q(x2)  >~6  (x^-x2),  has  been 

invoked  in  Section  1 1 -3 . Therefore,  one  might  suspect  that  < y(x^)y(x2)> 

would  be  fairly  sharp.  This  however,  turns  out  to  be  a false  speculation 

because  of  the  structural  rigidity,  so  that  the  spatial  correlation 

extends  over  the  entire  beam.  As  a typical  example,  let  us  compute  the 

cross-correlation  about  the  midpoint  of  beam.  Then,  with  the  use  of 
4 

I = 1,-,/n  , Equation  36  takes  the  following  form 
nn  1 1'  M 

< y (L/2  + Ax/2)y(L/2  - Ax/2)  > = [F-G]Kln,  (37) 


where 


F = 


odd 


I 

n = I 


COS 


( 


H7tAx  \ 
2L  / 


and 


G = 


even 

I 


n = 2 


The  square  bracket  of  Equation  37  has  the  peak  value  at  Ax  = 0,  and  falls 
off  to  zero  at  the  maximum  separation  Ax  = L.  It  is,  however,  more 
meaningful  to  examine  the  correlation  coefficient. 


V 


< y ( L/2  + Ax/2  ) y(L/2-Ax/2)> 


(L/2  + Ax/2 


)> 


(L/2  -Ax/2)  > 


I 


which  becomes 


V 


F - G 
F + G 


(38) 
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Since  F and  G are  always  positive,  it  follows  that  n < 1;  the  equality  is 
valid  only  when  Ax  = 0.  The  limiting  value  at  Ax  = L cannot  be  deduced 
directly  from  Equation  38,  for  there  it  is  indeterminate.  However,  it  is 
found  by  the  asymptotic  analysis  that  n = 1/2  at  Ax  = L.  The  presence 
of  spatial  correlation  over  the  entire  beam  is  depicted  in  Figure  4. 


In  the  remainder  of  this  section,  it  will  be  shown  that  the  prediction 
scheme  in  Equation  33  also  follows  from  Equation  34,  thus  demonstrating 
consistency  of  that  prediction.  To  the  first-term  truncation,  the 
variance  and  cross-correlation  are 


< 


> 


K I 


11  ’ 


< > (xi)  y ( *2)  > = 


(39) 


2 2 

where  < y (x-j ) > and  < y (x^)  > take  the  positive  and  negative  signs, 
respectively.  Now,  introduce  Equation  39  into  the  right-hand  side  of 
Equation  34.  After  averaging  over  the  beam,  one  obtains 

(VA«Ls=  (40) 


By  identifying  K 1^  with  2y2ms,  it  is  evident  that  (2  cos2(ttAx/2L)  - 1) 
plays  the  role  of  a beam-averaged  correlation  coefficient.  Further 
consolidation  of  Equation  40  by  factoring  out  K 1^  ( = 2y2rms)  yields 

(a’/a‘L,  ■ HH) • (41 , 

Since  2(1  - cos2  (ttAx/2L)  = 1 for  Ax  = L/2,  Equation  41  becomes 
identical  to  Equation  33,  thereby  establishing  the  internal  consistency 
of  Equation  33  by  way  of  an  alternate  derivation. 
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Ax 

L 


Figure  4.  Correlation  Coefficient 
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SECTION  VI 
CONCLUSIONS 

By  relating  the  linear  and  angular  responses  of  structural  vibration, 
qualitative  relationships  have  been  developed  to  predict  the  rms  amplitude 
of  angular  vibration  from  a prescribed  linear  vibration  response.  The 
prediction  schemes  have  been  proposed  for  the  true  angular  and  discrete 
angular  vibration  amplitudes.  Based  on  the  presently  available  test 
data,  the  overall  performance  of  predictions  is  within  +20%  of  the 
angular  vibration  measurement.  Although  crude,  such  a prediction 
capability  is  very  satisfying.  This  is  because  the  prediction  schemes 
are  intended  for  qualitative  purposes  and  the  present  test  data  are 
limited  only  to  gross  confirmation.  Since  the  present  investigation  is 
a first  attempt,  the  refinement  of  angular  vibration  predictions  will 
follow  as  more  accurate  vibration  data  become  available. 
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APPENDIX 

STOCHASTIC  ANALYSIS  OF  THE  TIMOSHENKO  BEAM 


For  a small  amplitude  vibration,  the  angle  of  rotation  may  be 

approximated  by  3y(x,t)/3x;  hence,  the  angular  velocity  and  angular 
2 3 2 

acceleration  are  3 y/3x3t  and  3 y/3x3t  , respectively.  The  Timoshenko 
beam  equation  which  takes  account  of  the  angular  damping  and  angular 
inertial  effects,  is  given  by 


a3y 

dx2dt 


J_  a4y 

0 dx2dy2 


f ( x,t) . 


(Al) 


The  last  two  terms  of  the  left-hand  side  modify  the  Bernoull i-Euler 
beam  Equation  1.  Here  s'  is  the  appropriate  damping  constant,  I is  the 
moment  of  inertia,  and  a is  the  cross-section  of  beam.  Although  the 
angular  inertial  term,  -(I/a)  3 y/3x  3t  , first  appears  in  Rayleigh's 
classical  beam  equation  (Reference  11),  it  is  Timoshenko  (Reference  1) 
who  has  investigated  both  the  angular  inertia  and  shearing  deformation 
effects  in  the  early  1920's;  hence  Equation  Al  is  referred  to  as  a 
Timoshenko  beam  equation.  More  recently,  Samuels  and  Eringen  (Reference 
12)  have  considered  the  most  general  Timoshenko  beam  equation  in 
search  for  a convergent  bending  stress  under  purely  random  excitations. 


For  the  simply  supported  case,  the  undamped  normal  modes  of  Equation 
Al  are  still  the  same  as  the  Bernoul 1 i -Euler  beam.  Therefore,  the 
stochastic  analysis  of  Section  II  can  be  repeated  verbatim  by  starting 
from  the  time-dependent  equations 

dn  yn  + 2/n  yn  + u>n  Yn  = 9n  9 (♦). 

2 2 

where  an  = 1 + (I/a)  ( nir/L ) and  Yn  = 8 + 8'  (n-rr/L)  (Compare  this  with 
the  Equation  following  Equations  2a  and  2b  in  the  main  text).  The 
linear  and  angular  variances  are  again  given  by  Equations  10  and  12, 
respectively,  except  for  the  modified  I ^ 

Ink  Fg  (tu)  Hp  (ct»)  H|<*  (oj)  dm,  (A2) 
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where  H'  (w)  = 2 * Under  the  white-noise  spectrum 

n + 2iy„  a)  - a a) 

n n n 

(Equation  15),  the  evaluation  of  Equation  A2  yields 

47r(Xn/qn>F0 

lnk"  (w2-wk2- 4yn/qn)2  + (4Xn/an)2(a;n2an-Xn2) 

which  for  n = k simplifies  to 

ttF0 

nn=  UiMWi?)  <«> 

Note  that  a does  not  appear  in  Equation  A3;  it  is  identical  to 
n o 

Equation  17  if  3 is  replaced  by  6 + 3'  ( nir/L ) . The  inclusion  of 
angular  effects  has  therefore  resulted  in  more  rapidly  decreasing  Inn 
with  respect  to  n,  I ~l/n**;  however,  it  has  no  effect  on  the 
prediction  schemes  derived  in  the  main  text. 
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